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Abstract 

We have extended the Interacting Gluon Model (IGM) to calculate diffractive mass spectra gen- 
erated in hadronic collisions. We show that it is possible to treat both diffractive and non-diffractive 
events on the same footing, in terms of gluon-gluon collisions. A systematic analysis of available data 
is performed. The energy dependence of diffractive mass spectra is addressed. They show a moderate 
narrowing at increasing energies. Predictions for LHC energies are presented. 
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1 Introduction 



In the last years, diffractive scattering processes have received increasing attention for several reasons. 
These processes may, for example, explain many features of particle production and, in particular, heavy 
flavour production jjj and Centauro events 0. They are also related to the large rapidity gap physics 
and the structure of the Pomeron . In a diffractive scattering, one of the incoming hadrons emerges 
from the collision only slightly deflected and there is a large rapidity gap between it and the other final 
state particles resulted from the other excited hadron. In some models diffraction is due to the Pomeron 
exchange but the exact nature of the Pomeron in QCD is not elucidated yet. The first test of a theory (or 
a model) of diffractive dissociation (DD) is the ability to properly describe the mass (Mx) distribution of 
diffractive systems, which has been measured in many experiments Q and parametrized as (M x )~ a with 
a ~ 1. Data presented in |I| were taken at the Tevatron collider (-^/s = 1.8 TeV). They allow us to make 
comparison with CERN || lower energies data and observe the energy dependence of the mass spectrum. 
In the old Regge theory, the assumption of Pomeron dominance implies that the mass spectrum behaves 
like l/M x and does not depend on the energy || whereas in ^ a slight deviation from this behaviour 
was reported. 



In this work we intend to study diffractive mass distributions using the Interacting Gluon Model 
(IGM) developed by us recently @, |§, [§• In particular, we are interested in the energy dependence of 
these distributions and their connection with inelasticity distributions Jl(j, || . One advantage of the IGM 
is that it was designed in such a way that the energy-momentum conservation is taken care of before 
all other dynamical aspects. This feature makes it very appropriate for the study of energy flow in high 
energy hadronic and nuclear reactions |?], |^, |^| and in cosmic ray studies ]To| . In particular, as shown in 
§,§, the IGM was very uscfull in analysing data and making predictions on the behaviour of inelasticities 
and leading particle spectra, including leading charm production Jll|] . The aim of the present work is to 
demonstrate that the main characteristics of the DD processes mentioned above emerge naturally from 
the standard IGM, enlarging profoundly its range of applications. 

In the following section we present for readers convenience the basic elements of the IGM. It is then 
applied to DD processes in Section 3 where we also carefully explain what DD means in terms of the 
IGM and how it differs from the conventional Regge Pomeron approach. Section 4 contains our numerical 
results and comparison with data and the last section is devoted to our conclusions. 



2 Interacting Gluon Model 

2.1 General ideas 

The Interacting Gluon Model is based on the following idea [[l2) : since about half of a hadron momentum 
is carried by gluons and since gluons interact more strongly than quarks, during a collision there is a 
separation of constituents. Valence quarks tend to be fast forming leading particles whereas gluons 
tend to be stopped in the central rapidity region. This picture is consistent with string formation and 
fragmentation as it is formulated in the Lund Model or in the Dual Parton Model. These models are 
based on the concept of string and were constructed to work at low average transverse momenta. When, 
in the late eighties, the energy of hadronic collisions increased by more than one order of magnitude it 
became necessary to incorporate the concept of parton and of hard and semi-hard collisions. The latter 
are collisions between partons at a moderate scale (Q 2 ~ (2GeV) 2 ) which however still allows for the 
use of perturbative QCD. The scattered partons form the so-called minijets. At yfs = 540 GeV the 
minijet cross section is already 25% of the total inelastic cross section. A new generation of models 
appeared using the concept of parton and the QCD parton model formalism and trying to understand 



minimum-bias multiparticle physics extending the parton model to semihard (moderate) energy scales. 
Among these models we mention those presented by Gaisser and Stanev Sjostrand p] , Wang |ll| 
and Geiger |16 . The IGM belongs to this class of models. In all these models one finds at a certain point 
expresssions of the type 

/ dxi I dx 2 f(xi,Q 2 ) f(x 2 ,Q 2 )cr(x 1 ,X2,Q 2 ) (1) 
JQ 2 J9l 

where f and a are parton momentum distribution and parton-parton elementary cross section respectively. 
Q 2 is the scale. Apart from some ambiguity in choosing the scale, these models have to face the problem 
that even at very high energies a significant part of a hadronic collision occurs at scales lower than the 
semihard one. At this point f and a are not very well known. The atitude taken in HIJING |l5|, in 
the Parton Cascade Model |16 and also in the IGM is to extrapolate these quantities to lower scales. 
These extrapolations can be continuously improved, especially in view of the advance of our knowledge 
on non-perturbative effects. There are, for example, models for distribution functions which work at 
scales as low as 0.3 GeV 2 JllJ. As for a one can compute non-perturbative effects in the context of 
an operator product expansion Inspite of these limitations these models have the advantage of 

dealing with partons and being thus prepaired to incorporate perturbative QCD in a natural way. This 
is welcome since perturbative processes are expected to be increasingly important at higher energies. 
Compared to the other models mentioned above, the IGM is simpler because it is designed to study 
energy flow and makes no attempt to calculate cross sections or to follow hadronization in great detail. 
This simplifies the calculations and avoids time consuming numerical simulations. The most important 
aspect of the IGM, shared with those models, is the assumption of multiple parton-parton incoherent 
scattering which is implicit in the Poissonian distribution of the number of parton-parton collisions (which 
is also used in refs. jli], [l3j [ll| ) used below. In going to lower resolution scales this independent collision 
approximation becomes questionable. On the other hand lattice QCD calculations jl9| in the strong 
coupling regime indicate that the typical correlation length of the soft gluon fields is around 0.2 fm. 
This is still smaller than the typical hadronic size. Therefore the independent collision approximation 
may still be a reasonable one. From the practical point of view, it was shown in that replacing the 
Poissonian distribution by a broader one does not affect the results significantly as long as some mass 
scale is introduced to cut off the very low x region. 



2.2 Formulation of the model 

The IGM is based on the assumed dominance of hadronic collisions by gluonic interactions [[l2| and can 
be summarized as follows (cf. Fig. la): 

(i) The two colliding hadrons arc represented by valence quarks carrying their quantum numbers 
(charges) plus the accompanying clouds of gluons (which represent also the sea qq pairs and therefore 
should be regarded as effective ones). 

(ii) In the course of a collision gluonic clouds interact strongly and form a gluonic central fireball (CF) 
located in the central region of the reaction. 

(Hi) The valence quarks (plus those gluons which did not interact) get excited and form leading jets 
(LJ's) (or beam jets) which then populate mainly the fragmentation regions of the reaction. 

It should be stressed that the IGM has been formulated originally in order to give the initial conditions 
for hydrodynamical models by providing in a dynamical way the so called inelasticity of the reaction 
understood usually as the (invariant) fraction of energy stored in the CF: 



K = y/x ■ y 



(2) 



with x and y being the fractions of the initial energy-momenta of the respective projectiles allocated 
to the CF. This variable and, in particular, its energy dependence is of vital importance in cosmic ray 
studies as it is the necessary ingredient allowing to deduce any elementary information from cosmic ray 
experiments p0| . 

According to the IGM, pairs of gluons collide and form "minifireballs" (MF's). In collisions at higher 
scales a minifireball is the same as a pair of minijets or jets. In the study of energy flow the details of 
fragmentation and hadron production are not important. Most of the MF's will be in the central region 
and we assume that they coalesce forming the CF. The probability to form a fireball carrying momentum 
fractions x and y of two colliding hadrons is defined as the sum over (an undefined number n of) MF's 

= 



x{x, y) = s ( x ~ X! niXi ) 5 ( y ~ X! n%Vi ) II - p (™ i ) 



(3) 



(all masses and transverse momenta are neglected in what follows). The distribution of the number of 
MF's is given by Pijii) for which we use Poisson distributions: 

Pin,) = ^exp{-n t ) (4) 

corresponding to independent production. Expressing the delta functions via Fourier integrals one can 
perform all summations, transform certain summations over Hi into integrals and make the replacement 
dfii/dxdy = ui(x, y) arriving at the general formula 

x{x,y) = X L— ■ (5) 
•exp (--i- [{y 2 )(x - (x)) 2 + (x 2 )(y - {y) f - 2(xy)(x - (x))(y (y))} } , 



where 



and 



D xy = {x 2 ){y 2 ) - {xyf 



(x n y m ) = fdxx n fdyy™^^), (6) 



io Jo 

with xo being a normalization factor defined by the condition that 

r2 



dx I dyx{x,y)Q(xy - K min ) = 1 (7) 
o Jo 

with K m i n = being the minimal inelasticity defined by the mass mo of the lightest possible CF. The, 
so called, spectral function oj(x,y) contains all the dynamical input of the IGM in the general form of 
(cf. |) 

u{Xi y) = G{x) G{y) Q ( xy _ K * m . n ) , (8) 

where G's denote the effective number of gluons from the corresponding projectiles (approximated by 
the respective gluonic structure functions) and a gg and a the gluonic and hadronic cross sections, respec- 
tively. These are basic formulae of the IGM from which all its applications are derived EG]. 



3 Diffractive Dissociation in the IGM approach 



It is now quite straightforward to extend the IGM also to diffractive dissociation processes. In Fig. lb we 
show schematically the IGM picture of a diffractive dissociation event. Unlike in Fig. la here only one of 




the protons looses fraction x of its original momentum and gets excited forming a LJ carrying ij, = 1 — x 
fraction of the initial momentum. The other one, which we shall call here the diffracted proton, looses 
only a fraction y of its momentum but otherwise remains intact [ p5| . In the standard IGM presented in 
the previous section we were computing the probability \{x,y) of depositing energy fractions x and y in 
the central region in the form of gluonic CF of mass M = ^Jxy s, which subsequently was decaying and 
producing particles. Here we shall be rather interested in the mass Mx, a new variable in our problem, 
which, as can be seen in Fig. lb, is just the invariant mass of a system composed of the CF and the 
leading jet formed by one of the colliding protons (we shall call it also diffractive mass). Denoting by El 
and Pl the energy and momentum of the upper (in Fig. lb ) proton and by W and P the energy and 
momentum of the CF, 

P L =E L = ^(1-x), P=^( x ~y), W=^-(x + y), (9) 
the energy Ex and momentum Px of the diffractive cluster are given by: 

E x =E L +W = + and P x = P L + P = ^ (1 - y) (10) 

leading to the following expressions for the mass of our diffractive cluster, Mx, and its rapidity, Yx'. 

(ii) 

where y/s is the invariant energy of the pp system. We are working in the cm frame of the incoming 
nucleons. All masses have been neglected p6[ . 

In the limit y — ► 1, the whole available energy is stored in Mx which remains then at rest, i.e., 
Yx = 0. For small values of y we have small masses Mx located at large rapidities Yx- In order to 
regard our process as being trully of the DD type we must assume that all gluons from the target proton 
participating in the collision (i.e., those emitted from the lower vertex in Fig. lb) have to form a colour 
singlet. Only then a large rapidity gap will form separating the diffracted proton (in the lower part 
of our Fig. 1) and the Mx system (in its upper part), which is the experimental requirement defining 
a diffractive event. Otherwise a colour string would develop, connecting the diffracted proton and the 
diffractive cluster, and would eventually decay, filling the rapidity gap with produced secondaries. In this 
way we are effectively introducing an object resembling closely (but by no means identical to) what is 
known as Pomeron (iP) and therefore in what follows we shall use this notion as a handy abbreviation 
of notation p7| . The probability of finding the "bunch" of gluons (that forms the Pomeron) mentioned 
above in a colour singlet configuration is not computable at a more fundamental level because the model is 
simple and does not involve colour quantum numbers explicitely. However this probability is represented 
indirectly in the cross section between the gluonic bunch and the other proton, the parameter a discussed 
below in eq. disj). 

In our approach the definition of the object IP is essentially only kinematical (3^], very much in the 
spirit of those used in all other works which deal with diffractive processes in the parton and/or string 
language [52], |3^, [34j. We shall therefore try to derive the whole M x dependence directly from the 
IGM. Our first step necessary to adapt the standard IGM to DD collisions will be the introduction a 
kinematical restriction in the formula for the moments eq. @ preventing gluons coming from the diffracted 
proton (and forming our object IP) to carry more energy than the one released in the diffractive system. 
Therefore we shall write it as 

(x n y m ) = fdxx n [ Vma *d y y m u(x,y), (13) 
Jo Jo 

where y ma x = — — and the meaning of the spectral function u>(x, y) remains the same as before. 



We can now calculate the diffractive mass distribution Mx using the x[ x iV) function by simply 
performing a change of variables (cf. eq.(|ll|)), 

dN f 1 ^ 



dM x 



dx / dy x(x, y) 6 (M x - sy) 9 (xy - K^ in ) 
o Jo 



M 



dx x [ x, — — 



(14) 



In the IGM j7| |j| the distribution x(x, y) is a wide gaussian in the variables x and y changing slowly with 
the energy y/s. Substituting now eq.(||) into eq.(|lj) we arrive at the following simple expression for the 
diffractive mass distribution: 

A AT 1 

,1 , , r ,w2 ., (15) 



dN 



where 



and 



H(M x ,s) 



dM% 



- ■ F(M x ,s)-H(M x ,s) 
s 



F(M x ,s) = exp 



Xa 



2tt/D 



■nj 



dx ■ 



2D,, 



M x I \ 
s 



(16) 



(17) 



• exp 



(y 2 )(x-(x)) 2 -2(xy)(x~(x)) 



M 



x 



-(y) 



The moments (q n ), q — x,y and n = 1,2 are given by ( p^| ) and are the only place where dynamical 
quantities like the gluonic and hadronic cross sections appear in the IGM. At this point we emphasize 
that we are all the time dealing with a proton- proton scattering. However, as was said above, we are 
in fact selecting a special class of events and therefore we must choose the correct dynamical inputs in 
the present situation, specially the gluon distribution inside the diffracted proton and the hadronic cross 
section a appearing in u. As a first approximation we shall take G F (y) = G p (y) = G(y) (cf. |p8| ), with 
G(x) = p(m + 1)^—^—, with m=5, the same already used by us before |J. The amount of the diffracted 
nucleon momentum, p, allocated specifically to the IP gluonic cluster and the hadronic cross section a are 
both unknown. However, they always appear as a ratio (£) of parameters in u> and different choices are 
possible. Just in order to make contact with the present knowledge about the Pomeron, we shall choose 



a(s) 



61n- 



(18) 



where so = 1 GeV 2 and a and b are parameters to be fixed from data analysis. As it will be seen, cr(s) 
turns out to be a very slowly varying function of y/s assuming values between 2.6 and 3.0 mb, which is 
a well accepted value for the Pomeron-proton cross section, and p ~ 0.05. 



Before performing a full numerical calculation let us estimate eqs.(|l5|) and (JT^) keeping only the most 
singular parts of the gluonic distributions used (i.e., G(x) ~ 1/x) and collecting all other factors in eq.(||) 
in a single parameter c. Let us first assume that the ratio of the cross sections does not depend 



on x and y. Neglecting all terms of the order of 
the moments calculated in eq. (|l3|) : 



and we arrive at the following expressions for 



2(x z 



In— £: 



(19) 



Ml 



Ml 



Notice that in all cases of interest (x ■ y) is much smaller than other moments (by a factor In 
least). It means that D xy ~ (x 2 )(y 2 ) and consequently 

( ^ - (v) ' 2 

F(M%,a) 



at 



exp 



exp 



2(y 2 



1 - c • In 



Mi 



c ■ In 



(22) 



and 



H(M 2 x ,s) 



Xo 



cfa; exp 



(*-(^)) 2 
2(a; 2 ) 



— #■ 



const • = const • 



leading to 



const 



dN 



M| • Jc-ln 



(23) 



- • H(M x ,s)-F(M x ,s) 



const 



exp 



c • In 



Mi 



1 - c • In 



Mi 



•In 



Mi 



(24) 



The expression above is governed by the jp- term. The other two terms have a weaker dependence on 
M x . They distort the main (jtt) curve in opposite directions and tend to compensate each other. It is 
therefore very interesting to note that even before choosing a very detailed form for the gluon distributions 
and hadronic cross sections we obtain analytically the typical shape of a diffractive spectrum. 



4 Comparison with experimental data 

In Fig. 2a we show our diffractive mass spectrum and compare it to experimental data from the CERN- 
ISR jj5), which are usually parametrized by the form jp-. These spectra were calculated with expression 
( |l5| ) with the same gluon distributions, cross sections and the mass parameter mo (mo = 350 MeV) used 
in previous works ||, As it can be seen, the agreement between our curves and data is reasonable. 
At large values of — experimental points start to flatten out, deviating from the jp- behaviour. This 
may be due to the contribution of non-diffractive events. We expect therefore some discrepancy between 
theory and experiment in this region. At very low values of and lower energies -y/s our model 
does not give a good description of data. Here again some discrepancy should be expected because we 
are neglecting all resonance effects. In principle a better agreement between theory and data could be 
achieved in this region. It would be enough to choose mo = 550 MeV keeping everything else as before. 
The result is shown in Fig. 2b. A discrepancy persists at very low masses at the lowest center of mass 
energy t/s. However this region is almost beyond the validity domain of the model. Instead of changing 
mo (which would also require a change in our previous works), we prefer to keep its usual value. Given the 
qualitative nature of the present work, we choose to be consistent with our previous works in prejudice 
of the quality of the fits. 

Fig. 3 and 4 show similar comparisons for y/s — 546 and 1800 GeV respectively. Data are from refs. 
and yj. Again we find reasonable agreement with experiment. In the low region at higher y/s 
the agreement is very good. All the curves above were obtained with a = 2.6 mb and b = 0.01 mb. 



As a straightforward extension of our calculation we now apply expression ( |15[) to the study of diffrac- 
tive pion-proton and kaon-proton scattering. We first consider the cases p+ir — > p+X and p+K — > p+X. 
This corresponds to replace the proton by a pion or a kaon in the upper line of Fig. lb, everything else 
remaining the same. We must also substitute the gluon distributions in the proton, G[x) 1 by the cor- 
responding gluon distributions in the pion and kaon, taken from ]36| ]. Here, for simplicity, we take 
G w (x) — G K (x). This is supported by an ACCMOR collaboration data analysis |3^| . We also assume 
that tr Pp = c P7r = cr PK . The comparison between our results and data from the EHS/NA22 collaboration 
]3l| is shown in Figs. 5a (pions) and 5b (kaons). We may also have diffracted mesons, which undergo 
reactions of the type n + p — > n + X and K + p — > K + X. This corresponds to replace the proton by a 
pion or a kaon in the lower line of Fig. lb, substituting also the corresponding gluon distributions. The 
comparison between our results and experimental data ]3S| is shown in Fig. 6a (pions) and 6b (kaons). 
As it can be seen, a good description of data is obtained. 

Let us consider now the energy dependence of our results. The CDF collaboration studied single 
diffractive events and found some energy dependence in the diffractive mass spectrum. This fact is 
illustrated by writing 

s d(J SD 1 / 9 kx 

vsd dM% (M|) 1+e 1 ] 

where the factor e which has been reported to be [^9| e = 0.121 ± 0.011 at y/s = 546 GeV and 
e = 0.103 ± 0.017 at y/s = 1800 GeV, respectively. Considering the error bars one might say that 
this value is just constant (and that there would be no indication of energy dependence), but a real (al- 
beit small one) change in e is not excluded. Therefore, if confirmed, it would mean that the distribution 
becomes slightly broader. 

In the IGM everything is from the beginning energy dependent and so should be the diffractive mass 
distribution. We start analysing the analytical approximation eqs. (|l5|,^2], |23|). In (S) we see that 
the s-dependence factorizes and the function H(M^-) has the same shape for all energies, the difference 
being only a multiplying factor (in numerical calculations this behaviour is slightly violated) . In ( |22] ) the 
s-dependence does not factorize and remains in the moments or, equivalently, in the variable c . F(Mj/) 
is a broad function with maximum value determined by the moment (y) which increases with the energy, 
making F to "rotate" in a way that it becomes higher at lower values of M\ and becomes deeper at larger 
values of M\. When we multiply H (which goes essntially like 1/Mj^) by F it becomes steeper. This 
behaviour of F and H is illustrated in Fig. 7a and 7b respectively, which show the numerical evaluation 
of eqs. (|l6l ) and (0). The result of the numerical evaluation of eq. (|l5| ) is presented in Fig. 8. There we 
show the energy dependence of our diffractive mass spectra in proton-proton scattering. Fig. 8a shows 
diffractive mass spectra for y/s — 23.5 GeV (solid lines), 44.6 GeV (dashed lines) and 62.4 GeV (dotted 
lines). Fig. 8b shows spectra at y/s = 0.54 TeV (solid lines), 0.9 TeV (dashed lines) and 1.8 TeV (dotted 
lines). Finally Fig. 8c shows our prediction for the diffractive mass spectrum at LHC (y/s = 14 TeV) 
compared to the Tevatron one. The spectra in Fig. 8a and 8b are the same of Figs. 2, 3 and 4. It 
is interesting to note that the behaviour that we find is not in contradiction with data. In all curves 
we observe a modest narrowing as the energy increases. This small effect means that the diffractive 
mass becomes a smaller fraction of the available energy y/s. In other words, the "diffractive inelasticity" 
decreases with energy and consequently the "diffracted leading particles" follow a harder xf spectrum. 
Physically, in the context of the IGM, this means that the deposited energy is increasing with y/s but it 
will be mostly released outside the phase space region that we are selecting. A measure of the "diffractive 
inelasticity" is the quantity £ = ^j L . Making a trivial change of variables in eq.(|l5|) we can calculate its 
average value (£) : 

(£>(*) = f^^i (26) 



where £ m in (= — ) and ^ ma x (= 0.1) are the same used in for the purpose of comparison. In 

Fig. 9 we plot (£) against y/s. As it can be seen (£) decreases with y/s not only because £ m i„ becomes 
smaller but also because ^ changes with the energy, falling faster. This qualitative behaviour of (£) 
is in agreement with the estimate of the same quantity extracted from cosmic ray data analysis p0| . 
Also shown in Fig. 9 is the quantity (£ £ ) (which has been discussed in [ [39"| in connection with the energy 
dependence of the single diffractive cross-section) for e = 0.08 (dashed lines) and e = 0.112 (dotted lines). 

The energy behaviour of is determined by the moments (|l^) and (||). Since cr gg (s) and G{x) 
are the same as in previous works, being thus fixed, the only source of uncertainty in the s-dependence 
of the results is in the ratio which is the only free parameter in the model. All curves presented 
above were obtained with the choice ( |l8| ) made for a and with p ~ 0.05. We have checked that for a 
stronger growth of cr Pp with y/s the energy behaviour of -jjjt might become even the opposite of the one 
found here, i.e., the diffractive mass distribution would become broader at higher energies. However these 
strongly s-dependent parametrizations do not give an acceptable description of the existing data and were 
therefore excluded. Considering what was said above one might think that we can discriminate between 
different Pomeron-proton cross sections and we could use this model to extract cr Pp from data. We stress 
however that, in this model, only the ratio ^ enters effectively in the calculations and it is impossi- 
ble to completely disentangle these two variables. In this sense, the almost constancy (with y/s) of cr Pp 
may be just an indication of some increase of p with y/s in a way that the ratio — "scales " with the energy. 



5 Conclusions 

To conclude, we have shown that: 

(i) The original IGM (|fl 0, 0) is reasonably successful in describing non-diffractive events. With 

JVf 2 

only two natural changes, namely the introduction of the kinematical cut-off y m ax = — — and the 
multiplication of the function wbya constant factor reflecting the essentially unknown combination 
of _ZP-proton cross section a and the amount of gluonic energy-momentum of the diffractive proton 
allocated to the object IP, p (in form of the ratio £), (with all other parameters kept as in previous 
applications of IGM) it turns out to be also able to provide a reasonable description of diffractive 
events [^lj and their energy dependence. We predict that, at higher energies, a narrowing of the 
M\ distribution may be observed. 

(it) As is obvious from Fig. lb, it provides a detailed description of the diffractive cluster (of mass 
Mx)i disentangling it in a natural way into gluonic cluster of mass M containing the majority 
of produced secondaries (which corresponds to CF in the original IGM) and leading jet (inside 
the diffractive system) with momentum fraction x l carrying quantum numbers of the diffractively 
excited projectile. This fact may be very uscfull for the cosmic ray applications of the IGM, in 
particular to studies using DD like those presented in j2| Q (where such disentanglement seems to 
be important and so far was introduced in an ad hoc way only) . 
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Figure Captions 



Fig. 1 IGM description of a proton-proton scattering : a) general case, b) with the formation of a diffractive 
system of invariant mass Afj. 

Fig. 2 a) Diffractive mass spectrum for pp collisions calculted with the IGM (eq.([l4|)) and compared with 
CERN-ISR data f|. b ) the same as a ) with m o = 550 MeV. 

Fig. 3 Diffractive mass spectrum for pp collisions calculted with the IGM (eq.([l4|)) and compared with 
CERN-SPS Collider data @. 

Fig. 4 Diffractive mass spectrum for pp collisions calculted with the IGM (eq.([l4|)) and compared with 
FERMILAB Tcvatron data g]. 

Fig. 5 a) Diffractive mass spectrum for p + n + — > p + X collisions calculted with the IGM (eq.(p^)) and 
compared with experimental data [38| . b) The same as a) for p + K + — > p + X collisions. 

Fig. 6 a) Diffractive mass spectrum for 7r + + p — ► n + + X collisions calculted with the IGM (eq.([l4]) and 
compared with experimental data |38). b) The same as a) for K + + p — > K + + X collisions. 

Fig. 7 a) Energy dependence of the function F(M|-,s) (|f6|). b) Energy dependence of the function 
H (M%, s) (0) calculted with the IGM. 

Fig. 8 a) Energy dependence of diffractive mass spectra calculated with the IGM (eq. (p^|)). The solid, 
dashed and dotted lines represent spectra at y/s — 23.5, 44.6 and 62.4 GeV respectively, b) The 
same as a) for y/s — 0.54 (solid lines), 0.90 (dashed lines) and 1.8 TeV (dotted lines), c) The same 
as a) for y/s = 1.8 (solid lines) and 14 TeV (dashed lines). 

Fig. 9 Energy dependence of the "diffractive inelasticity" (£) and of (£ £ ). 
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